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1.  Introduction 


Polarization  defines  the  phase  and  amplitude  relationships  between  the  various  components  of 
wave  motion  and  is  significant  in  all  technological  applications  based  on  wave  propagation,  such 
as  optics,  seismology,  telecommunications,  and  radar  science.  As  opposed  to  other  fields,  elastic 
wave  polarization  within  the  U.S.  Army  has  received  relatively  little  attention  due  to  the  general 
difficulty  in  evaluating  it  experimentally.  It  is,  however,  a  well-recognized  fact  that  the  ability  to 
measure  and  characterize  the  polarization  of  ultrasonic  elastic  waves  could  lead  to  the 
development  of  novel  structural  diagnostic  tools,  which  could  rely  on  the  sensitivity  of 
polarization  to  surface  roughness,  cracks,  temperature,  or  residual  stresses,  among  others. 

The  theoretical  study  of  the  polarization  of  Rayleigh  surface  waves  presented  in  Miklowitz  (1) 
suggests  its  estimation  as  an  alternative  method  for  surface  stress  estimation.  A  polarization 
parameter  can  be  defined  as  the  ratio  between  the  maximum  in-plane  and  out-of-plane 
displacement  components  and  can  be  directly  related  to  the  state  of  surface  stress  and  thus  be 
used  as  an  absolute  measurement  of  pre-stress.  Studies  on  the  acousto-elastic  effect  on  Rayleigh 
waves  in  a  homogeneous  material  include  the  work  of  Hirao  et  al.  (2)  and  Duquennoy  et  al.  ( 3 ). 
The  analysis  of  polarization  of  ultrasonic  guided  waves  is  even  more  limited  and  has  mostly 
focused  on  theoretical  aspects  related  to  the  description  of  the  wavefield  and  to  tuning  criteria  for 
the  excitation  of  specific  wave  modes  through  surface-mounted  transducers  (4-6).  Recent  work 
by  the  present  authors  has  shown  the  experimental  validation  of  the  polarized  components  using 
one-dimensional  scanning  laser  Doppler  vibrometry  (SLDY)  for  Rayleigh  waves  (7)  and 
extended  to  Lamb  waves  to  quantify  damage  in  metallic  plates  (8,  9). 

One  specific  component  that  has  direct  Army  significance  is  the  stiffener-plate  configuration,  or 
so-called  T-joint,  which  is  ubiquitous  in  aerospace  components  and  most  commonly  represented 
in  the  wing  airfoil-spar  and  bulkhead-stiffener  interfaces  found  in  Army  rotorcraft.  A  nominal 
fillet  is  generally  designed  at  the  interface  to  alleviate  potential  stress  concentrations,  which  can 
produce  initiatory  cracks  under  high-cycle  loading  conditions.  However,  these  T-joint 
configurations  provide  a  challenge  for  structural  health  monitoring  (SHM)  strategies  that  employ 
wave  propagation-based  interrogation  techniques  such  that  partial  transmission,  reflection,  and 
mode  conversion  are  modified  by  the  degree  of  the  fillet  radius,  excitation  frequency,  material 
constituency,  and  stiffener  cross-sectional  thickness  (10).  Since  the  success  of  this  proposed 
SHM  strategy  to  assess  damage  severity  largely  depends  on  the  quantification  of  the  mode 
conversion,  a  fundamental  understanding  of  the  wave  interaction  due  to  the  T-joint  presence  is 
necessary. 

Previous  work  has  been  conducted  on  understanding  T-joint  configurations  of  various  material 
constituencies  and  a  range  of  damage  quantification  techniques.  Modal  techniques  employed  by 
Ooijevaar  et  al.  (11)  and  Sundaresan  et  al.  (12)  used  scanning  laser  vibrometry  to  detect  damage 
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produced  in  composite  plate  with  two  T-shaped  stiffener  sections  and  induced  long  saw  cuts  on 
an  aluminum  wing  panel,  respectively.  In  terms  of  wave  propagation  techniques,  Greve  et  al. 
(13)  analyzed  Lamb  wave  interactions  of  a  cracked  welded  joint  and  rolled  and  plate  girders 
using  two-dimensional  (2-D)  and  three-dimensional  (3-D)  finite  element  modeling  (FEM). 
Simulations  and  experiments  showed  cracks  that  are  oriented  perpendicular  to  the  direction  of 
wave  propagation  cause  strong  reflections.  Oppenheim  et  al.  (10)  modeled  a  2-D  simulation  of 
Lamb  wave  transmission/reflection  at  a  T-joint  with  no  fillet  radius  using  Abaqus*  and  a  3-D 
model  of  a  plate  girder.  Experimental  measurements  using  piezoelectric  (PZT)  sensors  were 
located  on  the  web  in  pulse-echo  mode,  and  the  response  demonstrated  dominant  So  mode 
reflections.  Similarly,  Morvan  et  al.  (14)  measured  Lamb  waves  by  a  laser  vibrometer  in  a  T- 
joint  generated  by  a  contact  piezo-composite  transducer  and  normal  component  of  the  surface 
wave  displacement  of  the  plate  and  stiffener.  Kim  et  al.  (15)  analyzed  mode  conversion  of  T- 
joint  with  collocated  PZT  sensors  before  and  after  the  stiffener  using  a  2-D  simulation  under 
plane  strain  assumptions.  Other  studies  concluded,  with  experimental  evidence,  that  damage- 
sensitive  features  are  successfully  extracted  as  damage-induced  Lamb  wave  modes  have  unique 
characteristics  compared  to  those  modes  produced  by  the  stiffener  or  varying  thickness  (16, 17). 

The  initial  aim  of  this  report  is  to  investigate  whether  the  polarization  nature  of  Lamb  waves  may 
be  utilized  as  a  damage  localization  technique.  Specifically,  an  analytical  formulation  of  the 
Lamb  wave  elliptical  trajectory  orientation  is  developed.  Based  on  the  formulation,  two 
simplified  damage  localization  parameters  are  proposed,  and  finite  simulations  demonstrate  their 
potential  use.  The  techniques  are  then  applied  to  metallic  and  composite  plates  at  varying 
excitation  frequencies  with  different  geometries  and  damage  type,  and  a  3-D  SLDV  setup  is 
implemented  to  extract  the  orthogonal  displacement  components.  The  primary  advantage  of  the 
techniques  is  shown  to  reduce  the  required  signal  processing  while  enabling  high-resolution 
damage  localization.  A  secondary  objective  of  this  report  is  to  further  understand  the  mode 
conversion  mechanism  at  the  T-joint  interface  using  frequency- wave  number  filtering  techniques  per 
Ruzzene  (18),  both  as  a  function  of  the  fillet  radius  with  and  without  blind-side  damage  and  with 
experimental  validation  that  provides  high-fidelity  spatial  resolution  using  3-D  laser  vibrometry. 


2.  Analytical  Lamb  Wave  Polarization  Formulation 


The  particle  trajectory  at  a  given  location  is  calculated  based  on  the  analytical  Lamb  wave  model 
provided  by  Raghavan  et  al.  (19).  For  the  given  formulation,  the  displacement  is  initially 
described  as  a  superposition  of  the  fundamental  symmetric  (S)  and  anti-symmetric  (A)  Lamb 
wave  modes  at  a  frequency  lower  than  the  first  anti-symmetric  mode  cutoff.  The  far  field 
approximations  are  used  and  written  in  the  following  compact  form  (1): 


*  ABAQUS  is  a  registered  trademark  of  ABAQUS,  Inc. 
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and 


where 


u^r,  co)  =  Afe'^s  n/T>  +  Af  e'^A  ”/2) 
u3(r,  co)  =  A|e1(ps  +  A3e1<pA  , 


<Ps  =  ~ksr  +  -  ,  <pA  =  ~kAr  +  ~ , 


(1) 

(2) 

(3) 


such  that  uj  and  113  are  the  in-plane  and  out-of-plane  displacements,  r  is  the  radial  distance  away 
from  the  piezoelectric  source,  to  is  the  frequency,  A  is  the  amplitude  of  the  respective  displacement 
component  and  mode,  and  k  is  the  wave  number.  Normalization  of  the  displacements  by 


Real(uj) 

u>=~m~  0  =  U) 


results  in 


and 


-  _  A{  sin((p5)+A^  sinC^) 

Ui  —  1 - 

(Af)2+(4)2+2Af4  cos(<p) 


-  _  Af  sin(<p5)+A^sin(^) 

u3  ~  1 - 

(A§)2+(Af)2+2A§Af  cos  (<p) 


Using  trigonometric  identities  within  equations  5  and  6  leads  to  the  following: 

^  (Af  sin  (ps  +  Af  sin  cps  cos  <p  +  Af  cos  <ps  sin  cp)  , 


and 


where 


u3  =  -  (A§  cos  cps  +  A3  cos  cps  cos  cp  —  A3  sin  cps  sin  cp)  , 


T  =  V  (Af)2  +  (Af)2  +  2AfA3  cos  <p  . 
In  matrix  form,  equations  7  and  8  are  rewritten  as  follows: 


(Ui]=rBi1  B12lfsin<p5| 
(u3j  [b21  B22J  (cos  (ps)  ’ 


where 


and 


Bn  =  p  (Af  +  Af  cos  cp),  B12  =  p  Af  sin  <p 


-Laa, 


B2i  =  “  p  Af  sin  (p  ,  B22  =  -  (A f  +  Af  cos  cp)  . 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 
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Solving  the  matrix  in  equation  10  results  in  the  following: 


sin  (ps=\  (uiB22  ~  u3B12 )  ,  (12) 

and 

cos<ps  =  ^(-u1B21+u3B11'),  (13) 

where  B  =  B11B22  —  B12B21.  Squaring  and  summing  equations  12  and  13  produces 

—  [(Ui522  —  u3B12)2  +  (— iti£?21  +  u3  51;l)2]  =  1  .  (14) 

Upon  assigning  the  following  parameters, 

a  =  j  =  =,  (15) 

JB22+B12 

b  =  TJ=,  (16) 

Jb?i+b|i 

and 

cos  D  =  ,  =  (17) 

JB22+B12 

Applying  equations  15-17  to  equation  14  results  in  a  rotated  elliptical  equation  (figure  1)  as 
follows: 


COSO  _ 

(«i  — - u3 


sinQ^n  ,  sin  £2 
— )  +  (“I  -Z- 


+  U3 


cos  £2. 


)2  =  i  . 


b 


(18) 


Physically,  it  can  be  shown  that  the  representation  of  the  polarized  components  of  the  guided 
waves,  u-l  and  u3,  yield  an  elliptical,  rotated  profile  of  the  particle  trajectory. 


Figure  1.  Lamb  wave  elliptical  trajectory 

profile  with  identified  parameters. 
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Figure  2a  illustrates  the  Lamb  wave  polarization  ellipses  as  a  function  of  r,  the  radial  distance 
away  from  the  piezoelectric  source.  Figure  2b  shows  the  projections  of  these  ellipses  in  (ut,  r) 
plane. 


Figure  2.  Image  showing  (a)  normalized  elliptical  trajectory  as  function  of  the  distance  away  from  the  source  and 
(b)  planform  view  of  the  in-plane  polarized  component  as  function  of  the  distance  away  from  the  source. 

2.1  Damage  Localization 

Two  damage  localization  parameters  based  on  the  elliptical  trajectory  orientation  (ETO)  are 
proposed.  The  first  is  based  on  the  observation  that  the  period  of  the  orientation  (when  shown  as 
a  function  of  the  distance  of  the  piezoelectric  actuator)  varies  in  relationship  to  the  period  of  the 
incident  wave.  The  localization  parameter,  denoted  as  IIeto,  may  be  formulated  as  a  ratio  of  the 
resulting  period,  Qj,  and  the  incident  period  such  that: 

Tn.  ( Ar,co ) 

=  <19> 

where  T  is  the  period  of  the  elliptical  orientation  as  function  of  r  and  subscripts  i  and  j  are 
incident  and  resulting  periods,  respectively.  The  parameter  is  a  function  of  Ar  as  a  result  of  the 
period  calculated  from  the  orientation  changing  over  the  radial.  This  formulation  has  the  distinct 
advantage  of  relying  on  time  domain  spatial-temporal  response  and  thereby  eliminates  the 
domain  signal  processing  that  previous  damage  localization  parameters  have  commonly 
employed  (20). 

A  second  proposed  damage  indicator,  denoted  for  convenience  as  TETo,  is  based  on  the 
amplitude  differential  of  the  periodic  elliptical  orientation,  5;,  as  a  function  of  the  distance  away 
from  the  source  such  that 


and 


5j  (Ar,  co)  =  fli+1  —  fl j 


(20) 
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^ET0(Ar,a))  , 


(21) 


where  the  'Feto  parameter  is  normalized  by  the  maximum  periodic  amplitude  differential. 

2.2  Damage  Quantification 

This  analysis  is  restricted  to  the  frequency  range  of  the  fundamental  symmetric  and 
antisymmetric  Lamb  waves,  So  and  Ao,  respectively.  As  illustrated  in  figure  3,  it  is  assumed  that 
a  pure  So  or  A0  mode  is  generated  to  interrogate  the  structure.  The  interaction  of  an  incident  So 
mode  with  a  vertical  stiffener  causes  a  partial  So  reflection  and  partial  mode  conversion  into  the 
Ao  mode  while  simultaneously  producing  a  partial  So  and  Ao  mode  transmission  into  the  stiffener 
and  the  plate.  As  previously  described  by  Ayers  et  al.  (21),  upon  filtering  the  data  in  the 
frequency-wave  number  domain  and  windowing  the  data  in  the  time-space  domain,  w(x,  co)  is 
the  filtered  response  in  the  spatial-frequency  domain  of  a  single  extracted  mode  at  frequency  to. 


Laser  vibrometry  techniques  allow  for  various  points  to  be  scanned  over  the  distance,  L.  From 
the  measured  responses,  the  value  of  the  mode  coefficients  may  be  estimated  over  a  spatial 
domain.  This  allows  the  estimation  of  its  spatial-independence  and  more  importantly  of  a 
spatially  averaged  value,  which  is  less  affected  by  amplitude  fluctuations  due  to  noise  or  by 
general  trends  associated  with  geometrical  spreading.  In  compact  matrix  notation,  the  mode 
coefficients  may  be  expressed  by  the  following  for  an  interrogating  So  mode  but  are  likewise 
easily  applied  for  an  interrogating  Ao  mode  as  follows: 


n„(w)  =  - 


c o)\dx 

^ L 

*2 


(22) 
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where 


n,(w)  = 


R“ 
R w 

J^UW 

S o 


rriU 

1S 0 


rjiUW 

1S 0 


C“ 

Ao 

c; 

Ao 

(JUW 

So 


(23) 


%(w)  = 


r(r) 


(0 


rco 


(0 


w(r) 


ws0  Ws0 

Pl’f+wy  p£! ?+(*%)'- 
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and 


0((w)  = 


u 


(0 


W' 


,(0 


(24) 


(25) 


Q,(x)  = 


x1  x3  x1 


(26) 


where  R,  T,  and  C  parameters  represent  the  reflection,  transmission,  and  converted  coefficients,  u 
and  w  are  the  in-plane  and  out-of-plane  displacement  components,  respectively,  and  x,  values  are 
the  integration  spatial  coordinates.  It  is  important  to  note  that  this  spatial  integration  method 
departs  from  existing  literature,  where  the  mode  amplitudes  are  taken  at  a  fixed  location  (21). 


3.  Numerical  Results 


3.1  Damage  Localization 

A  simplified  plane  strain  finite  element  model  is  employed  to  evaluate  the  feasibility  of  the 
damage  localization  parameters.  Specifically,  an  isotropic,  homogeneous  plate  with  varying 
notch  types  at  the  mid-span  is  developed  using  the  commercial  software  Abaqus.  The  plate 
dimensions  are  2  x  1CT3  m  in  height,  5  x  10_1  m  in  width,  and  1.0  m  in  length.  The  material 
properties  reflect  those  of  nominal  A1  7075-T6,  where  the  density  p  =  2800  kg/m  ,  elastic 
modulus  E  =  72  GPa,  and  Poisson’s  ratio  v  =  0.33.  A  longitudinal  uniform  pressure  tip  load  is 
applied  at  one  end  of  the  plate,  whereas  the  opposite  end  is  fixed.  The  symmetry  of  the  applied 
in-plane  load  ensures  that  only  the  So  mode  is  excited,  which  is  later  verified  in  the  frequency- 
wave  number  domain  during  data  post-processing.  The  equations  of  motion  are  solved  by  the 
explicit  central-difference  time  integration  rule,  which  is  suitable  for  large  models  with  short 
dynamic  response.  The  central-difference  operator  is  conditionally  stable,  and  in  order  to  have 
an  accurate  solution,  the  integration  time  step  At  is  chosen  such  that  At  =  fo/20,  where  fo  is  the 
highest  frequency  of  interest,  measured  in  Hz  (22). 
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For  a  given  excitation  of  200  kHz  in  figure  4a,  the  incident  period  over  the  radial  distance  from 
0  to  0.2  m  is  0.013  m  (figure  4b),  which  includes  the  incident  and  reflected  fundamental 
symmetric  modes.  In  contrast,  from  figure  4b,  the  resulting  period  over  the  (0.2-0.5-m)  radial 
distance  is  0.06  m,  which  is  the  summation  of  the  incident  and  reflected  fundamental  symmetric 
mode,  So,  and  an  additional  converted  antisymmetric  mode,  Ao.  Applying  equation  19,  the 
consequent  Heto  is  shown  in  figure  5,  where  the  amplitude  is  unity  except  at  the  periodic 
discontinuities  at  r  =  0.2,  0.5,  and  0.8  m.  These  discrete  periodic  discontinuities  indicate  the 
presence  of  mode  conversion  and  consequent  damage  localization  without  additional  domain 
signal  processing  algorithms. 
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Figure  4.  Image  showing  (a)  in-plane  spatial-temporal  domain  of  surface  response  with  indicated  notch  damage 
location  at  r  =  0.5  m  and  (b)  elliptical  trajectory  orientation  as  a  function  of  radial  distance. 


Figure  5.  Damage  localization  parameter,  IIeto,  as  a  function  of 
radial  distance. 
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3.2  Damage  Quantification 

The  mode  coefficient  formulation  is  applied  to  numerical  parameterization  for  the  undamaged 
and  damaged  configurations  using  simplified  2-D  plane  strain,  where  mode  conversion 
estimation  is  shown  as  a  function  of  the  fillet  radius,  excitation  frequency,  and  stiffener  cross- 
sectional  thickness  (table  1  and  figure  6). 

Table  1 .  Test  matrix  of  2-D  FEM  parameterization. 


Set 

Signal 

Plate-Stiffener  Geometry 

Notch  Damage  Geometry 

Frequency 

(kHz) 

Fillet  Radius 
(mm) 

Stiffener 

Thickness 

(mm) 

Depth 

(mm) 

Orientation 

o 

1 

200.0 

0.0:3.0:12.0 

3.0 

— 

— 

2 

200.0 

0.0:3.0:12.0 

3.0 

3.0 

0.0 

3 

200.0 

6.0 

3.0 

0.0:2. 0:8.0 

0.0 

4 

200.0 

6.0 

3.0 

3.0 

0.0:15.0:75.0 

(a) _ (b) 


Figure  6.  Schematic  of  2-D  FEM  T-joint  showing  the  following:  (a)  profile  with  identified  boundary  conditions 
and  (b)  magnified  view  of  FEM  parameters  at  fillet  with  notch  damage. 

From  figures  7a  and  b,  the  mode  coefficients  from  a  no-damage  configuration  vary  parabolically 
as  a  function  of  the  fillet  radius.  With  a  constant  damage  depth  in  figure  8a,  the  RSo  mode 
coefficients  varies  linearly  as  a  function  of  the  fillet  radius.  As  the  fillet  radius  increases  (>6  mm), 
the  mode  coefficients  from  undamaged  to  damaged  are  approximately  equivalent  (figure  8b).  The 
thickness  of  the  fillet  radius  (>6  mm)  requires  that  the  damage  depth  does  not  penetrate  the  plate 
thickness,  and  sensitivity  is  negligible.  The  RSo  mode  coefficient  is  sensitive  to  the  increasing 
damage  depth,  whereas  the  Cao  mode  coefficient  shows  perturbation  (>4  mm)  once  the  damage 
depth  penetrates  the  plate  thickness  (figures  9a  and  b).  The  RSo  and  ^Ao  mode  coefficients 
demonstrate  a  sinusoidal  profile  with  the  notch  orientation  (figures  10a  and  b). 
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Figure  7.  Variation  of  fillet  radius  without  blind  damage  showing  (a)  RSo  mode  coefficient  and  (b)  Cao 
mode  coefficient. 


Figure  8.  Variation  of  fillet  radius  with  blind  damage  of  constant  depth  and  orientation  showing  (a) 
RSo  mode  coefficient  and  (b)  Cao  mode  coefficient. 


Figure  9.  Variation  of  blind  damage  depth  with  constant  fillet  radius  showing  (a)  RSo  mode  coefficient 
and  (b)  CAo  mode  coefficient. 
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Figure  10.  Variation  of  blind  damage  orientation  with  constant  fillet  radius  showing  (a)  RSo  mode 
coefficient  and  (b)  CAo  mode  coefficient. 


4.  Experimental  Results 


4.1  Damage  Localization 

The  damage  localization  parameters  were  experimentally  verified  by  testing  homogenous 
metallic  and  heterogeneous  composite  specimens  with  varying  types  of  induced  damage.  For 
this  report,  one  of  the  experimental  results  is  presented,  specifically  dynamic  excitation  of  a  low- 
cycle,  fatigue-induced  crack  on  a  3.175-mm-thick  6061-T6  aluminum  plate  cut  into  a  dog-bone 
shape  (figure  11a).  A  single  6.35-mm-diameter  x  0.254-mm-thick  PZT  transducer  was  bonded 
approximately  on  the  longitudinal  centerline  of  plate  from  the  hole  (center  of  hole  to  center  of 
PZT)  using  M-Bond  200  adhesive.  For  the  experimental  testing,  a  Hamming-windowed, 

5  1/2-cycle  sine  burst  excitation  signal  was  generated  using  an  Agilent  33120A  arbitrary 
waveform  generator.  This  signal  was  amplified  using  a  Krone-Hite  7500  wideband  power 
amplifier.  Data  collection  was  performed  using  a  Polytec  PSV-400-3D-M  SLDV  system 
designed  for  full-field  vibration  measurements  at  frequencies  up  to  1  MHz  (figure  lib).  A  dense 
grid  of  200  vertical  data  points  was  applied  for  measuring  the  wavefield  response,  with  an 
average  resolution  of  0.7  mm. 

For  a  given  tone  burst  of  100  kHz,  the  resulting  elliptical  orientation  is  presented  in  figure  12a. 
The  incident  period  over  the  radial  distance  from  0.02  to  0.04  m  consists  of  incident  and 
reflected  So  modes,  with  a  converted  Ao.  The  resulting  period  over  the  0.02-0.042-m  distance  is 
the  summation  of  the  incident  and  reflected  fundamental  symmetric  mode  and  an  additional 
transmitted  Ao  mode.  Applying  equation  19,  the  consequent  nETO  is  shown  in  figure  12b.  The 
approximate  rise  in  slope  and  effective  discontinuity  at  r  =  0.042  m  is  a  clear  indication  of  the 
mode  conversion  and  identifies  the  crack  presence. 
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Figure  11.  Images  showing  (a)  fatigued  aluminum  dog  bone  specimen  with  53-mm  crack  and  (b)  3-D  SLDV  test 
setup  with  data  acquisition  (DAQ)  and  test  stand. 


Figure  12.  Images  showing  (a)  elliptical  orientation  of  aluminum  dog-bone  specimen  and  (b)  damage  localization 
parameter,  nET0,  as  a  function  of  radial  distance. 


4.2  Damage  Quantification 

The  trends  from  the  FE  parameterization  are  investigated  by  using  a  3-D  scanning  laser 
vibrometry  test  setup  (figure  13a),  which  consists  of  an  aluminum  T-joint  for  selected  fillet 
radius  (0.6  cm),  excitation  frequency  (100  kHz)  and  stiffener  thickness  (0.3  cm).  The  setup 
enables  3-D  surface  measurement  of  the  front  and  back  sides  of  the  plate,  plate-stiffener  radius, 
and  stiffener  (figure  13b). 
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Figure  13.  Experimental  validation  using  Polytec  3-D  vibrometer  showing  (a)  test  setup  and 
(b)  aluminum  T-joint  specimen. 


To  ensure  Lamb  wave  phenomena,  figure  14a  illustrates  the  frequency  peaks  of  the  incident  Ao 
mode  of  the  front  face,  which  aligns  with  the  analytical  Ao  mode,  while  figure  14b  shows  the 
back  face  frequency  content,  where  mode  conversion  occurs  and  the  So  mode  is  introduced. 
Figure  15a  demonstrates  the  full  wave  field  of  the  filtered  out-of-plane  Ao  mode  along  the  front 
face,  and  figure  15b  depicts  the  in-plane  converted  So  mode.  The  magnitude  of  the  fillet  radius 
in  figure  15b  is  three  orders  below  that  of  the  front  and  stiffener  surfaces.  Using  data  whose 
spatial  points  originate  at  Y-axis  0.14  m  in  figure  15a  and  runs  along  the  front,  back,  and 
stiffener  faces,  table  2  contains  the  reflection,  transmission,  and  conversion  coefficients. 
Particular  attention  is  given  to  the  difference  in  order  of  magnitude  between  the  stiffener  and 
front  and  back  faces,  indicating  that  a  significant  amount  of  energy  is  propagating  into  the 
stiffener.  The  relative  order  of  magnitude,  -0.02-0. 06  along  the  front  face  for  the  reflection  and 
conversion  mode  coefficients,  corroborates  with  the  results  from  the  undamaged  FE  trade  studies 
shown  in  figure  3. 
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Figure  14.  Incident  in-plane  frequency  content  showing  (a)  front  face  and  (b)  back  face.  Excitation 
frequency  of  modulated  tone  burst  is  100  kHz. 


Figure  15.  Time  history  snapshot  at  87.8  ps  along  front  face-stiffener  surfaces  showing  (a)  filtered 
incident  out-of-plane  A0  mode  and  (b)  filtered  converted  in-plane  S0  mode. 


Table  2.  Integrated  mode  coefficients  along  center  line  of  three  surfaces. 


Section 

Reflection 

Coefficients 

Transmission 

Coefficients 

Conversion 

Coefficients 

Raou 

n  w 

*V\0 

n  u  w 

K\0 

rp  U 

aA0 

rp  W 

AA0 

rp  UW 

1 A0 

Csou 

/-I  w 
<-S0 

p  UW 

<-S0 

Front  face 

0.054 

0.059 

0.058 

— 

— 

— 

0.053 

0.024 

0.039 

Back  face 

— 

— 

— 

0.012 

0.019 

0.017 

0.014 

0.002 

0.009 

Stiffener 

— 

— 

— 

0.193 

0.203 

0.199 

0.118 

0.019 

0.074 
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5.  Conclusions 


An  analytical  formulation  of  the  elliptical  trajectory  is  provided,  and  the  technique  using  the 
orientation  of  the  elliptical  trajectory  as  a  damage  localization  parameter  is  investigated  on  an 
aluminum  specimen  with  an  induced  low-cycle  fatigue  crack.  The  crack  presence  was 
successfully  localized,  and  future  work  will  consist  of  extending  the  technique  to  full  damage 
characterization,  which  will  include  damage  quantification  and  complete  topology  of  the 
damage.  In  addition,  results  not  presented  herein  will  investigate  the  damage  localization 
techniques  to  more  complex  test  specimens  that  include  geometrical  complexity  and  material 
anisotropy. 

In  terms  of  damage  quantification,  as  the  fillet  radius  increases,  the  mode  coefficients  from 
undamaged  to  damaged  are  approximately  equivalent.  The  thickness  of  the  fillet  radius  (>6  mm) 
does  not  allow  the  damage  depth  to  penetrate  the  plate  thickness,  and  sensitivity  to  the  relative 
damage  is  therefore  negligible.  In  addition,  the  RSo  and  Cao  mode  coefficients  demonstrate  a 
sinusoidal  profile  with  the  notch  orientation  such  that  determining  the  orientation  based  on  the 
magnitude  of  the  mode  coefficient  is  limited.  Experimentally,  the  3-D  laser  vibrometry  results 
indicate  that  the  transmission  coefficients  are  an  order  of  magnitude  greater  within  the  stiffener 
section  than  the  front  and  back  faces. 
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